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On using Randomized Hough Transform for Detecting Radio Systems with Periodic 

Emission Pattern 



As indicated in earlier sections, when radio networks encounter other devices that emit 
energy (and therefore use shared radio resources) in their vicinity, it is desirable to 
characterize the radio resource usage patterns of these other devices. Such a 
characterization of the usage patterns results in the identification of opportunities for the 
radio networks. Other devices referred to previously includes radars, which are primary 
emitters, or other radio networks, which are secondary emitters. 

In this section we will examine the use of Hough Transform for the detection of radar 
pulses as an example for any type of radio signals that create periodic patterns. We will 
use a version of the Hough Transform, known as Randomized Hough Transform (RHT) 
to detect the parameters of helixes wrapped around cylinders, as explained later in the 
section. The Hough Transform [] has been studied in image processing literature for 
detection of patterns such as lines, circles and ellipses in binary images. The effectiveness 
of Hough Transform in detecting patterns in "imporfoot" data with many overlaying 
patterns and random noise is proven []. In the presence of outliers, the Hough Transform 
is more robust than least squares estimation []. 

In brief, the Hough Transform is used to transform data from image space to an 
accumulator (or histogram) in parameter space, as shown in the Figure 1 below. 
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(a) 0>) 
Figure 1 : Hough Transform used to detect straight 
lines (a) image space and (b) parameter space. 
The image space is represented by (x, y), whereas, the parameter space is represented by 
(slope, intercept), that is (m, c). For each point in the image space (e.g. p and q), a line is 
generated in the parameter space as shown. The parameter space can be seen as a two 
dimensional histogram. A peak (r) in the parameter space corresponds to a line in the 
image space. The Hough Transform is robust because in the image space, a collection of 
collinear points is enough to result in a peak in the parameter space. However, it has the 
drawback that the parameter space could require large amount of memory in the 
computer. To address this drawback the RHT was developed. The RHT as applied to 
straight-line detection, results in randomly picking pairs of points and computing and 



accumulating a parameter (for instance, slope). When enough confidence in the peak is 
achieved, the process stops, thus reducing both memory and processing time. 



The use of Hough Transform for radar pulse detection was first studied in []. The original 
radar pulse train is a 1-D signal. The authors have used 1-D to 2-D transformation (like a 
raster scan) and then applied the Hough Transform to detect straight lines, which 
correspond to pulse trains. Furthermore, they have computed the noise floor. We extend 
their work by first transforming the 1-D signal to a 3-D helical signal, and apply RHT to 
it. 

A helix may be represented by the following parametric equations: 

7(0 = cos(fi)0 (1) 

Z (<) = < 

This helix is cylindrical (as opposed to the more general elliptical) and has unit radius. 




Figure 2: A helix given by the equation 1 above (for ©=1) 
Based on the parameter a> a new helix can be generated that wraps around the cylinder 
more slowly as <o decreases. In Figure 2, the points (red *) on the blue helix themselves 
form a helix, with an <o value less than one. Given two points on the helix Po (*o, Yo, ^) 
and Pi (x,, y,, z,), the parameter a can be given by the following equation: 



atan(-^-)-otan(-^-) 
X\ Xq 

Z\"Zo 



.(2) 



If the two points P 0 and Pi are inside one twirl of the helix, then ©works out to be 1. The 
length of the line segment given by one twirl of helix is 



7 = 2***. 



1 + 



f-1 



(3) 



Results 



Let us represent the location (time-of-arrival) of the radar pulse train with the vector L p . 
ForL p / = [9, 59, 109, 159, 209, 259, 309, 359], the following go histogram is obtained. 




Figure 3: Histogram of o for L pl 



For this case, © = 0.1 16. Now let us consider the case where there are two pulse trains 
that are multiplexed and represented by L p2 = [9, 20, 59, 60, 100, 109, 140, 159, 180, 209, 
220, 259, 260, 300, 309, 340, 359]. The following histogram is obtained. Note that co = 1 
corresponds to points on the helix within one twirl. 
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Figure 4: Histogram of© for L p2 . 
Potential Extension 

The method described in this document may be associated with the alternative way of 
detection periodic interferences, which is based on the autocorrelation of measurements 
in time, We expect that both alternatives show advantages and disadvantages in different 
scenarios, and a combination of both may therefore result in the most precise 
identification of other radio systems. 
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